The role of Coulomb interaction in fragmentation 
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We examine the impact of adding a Coulomb term to a constrained system of 147 particles inter- 
acting via a Lennard- Jones (I/J) potential, finding that the inclusion of the coulombic interaction 
produces a shift, but no qualitative changes in the thermodynamical properties of the system. We 
also performed the systems characterization from a morphological point of view. 
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The analysis of the caloric curves (CC) of small sys- 
tems is attracting the attention of physicists in different 
areas. One of the most active ones is precisely nuclear 
physics because of the possibility of the occurrence of 
phase transitions in multifragmentation experiments. A 
wealth of theoretical and experimental work has flour- 
ished in recent years with, in some cases contradictory 
results. In particular a debate has recently arisen re- 
garding the effect of Coulomb interaction on the proper- 
ties of the CC. Whereas some authors references [l|, Q 
claim that there exist a loop in the CC denoting a neg- 
ative thermal response function for systems as large as 
200 particles it has also been claimed that an upper limit 
in the mass exist {A = 60) 0, 0| above which no such a 
loop can exist. In Ref.0 it is stated that the presence 
of Coulomb forbids the phase transition for systems as 
large as A ^ 200, Z = 82 (a proton fraction of 0.4). On 
the other hand in Ref.Q they propose a model in which 

is never negative. 

In this communication we will show that, when dealing 
with finite constrained systems interacting via a Lennard- 
Jones 6 — 12 potential which has been shown to display 
negative values of the specific heat [or more generally 
speaking, the thermal response function (TRF)] for low 
values of the density p], the addition of the Coulomb 
interaction term preserve the shape of the caloric curve 
and gives rise to a displacement of the location of the loop 
together with its flattening. These results are consistent 
with lattice models calculations including Coulomb in- 
teraction Q. It is worth to mention that the relevance 
of the analysis of LJ systems relies on the fact that the 
EOS of the LJ and the one assumed to describe nuclear 
matter are quite similar [9lll0|. 

The system under study is composed by a gas of 147 
particles confined in a spherical box, defined by the 



Hamiltonian H = K + Vlj + Vc 



coulomb 



where 



K is the kinetic energy and the short-range term of the 
interaction potential is given by Vlj = Sw(ry) with 

v{r,j) = 4e[;gy'' - ;^^] if (0 < T < r,) and other- 



wise. In this work we took a cut-off radius Vc — 3. Oct, 
and adopted adimensional units for the energy, length 
and time such that e — a — 1, to — y'(r^m/48e. 

The role of Vwaiis is to constrain the particles inside 
an spherical container. The considered external potential 
behaves like Vwaii ~ (r — rwaii)~^^ with a cut-off distance 
fcut = Ic, where it smoothly became zero along with its 
first derivative. 



The addition of the coulombic interaction Vr 



coulomb 1^ 
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the hamiltonian has been done with two different ap- 
proaches that will be henceforth referenced as the homo- 
geneous and the inhomogeneous one. 

The homogeneous case, which could be useful in the 
study of metallic clusters and has already been intro- 
duced for the study of nuclear fragmentation jQj, the 
coulombic force is included between all particles of the 
fluid, Vcoukmib = Q^l^ij, with = 0.055, wich was ob- 
tained in [9j by comparing the binding energy formula 
for argon balls with the nuclear drop mass formula as- 
suming the proton fraction to be 0.4 (roughly speaking, 
it has been adjusted so that liquid argon drops having 

300 atoms are unstable under fission). 

In the second case, i.e. the inhomogeneous case, the 
Coulomb term is only present between a subset {Z) of 
the {A) particles which will be renamed as "protons" (the 
other A — Z will therefore called " neutrons" ) . Taking the 
same proton fraction as in the previous case (^ = 0.4), 
we adjusted the strength of the Coulomb term so that we 
have the same total energy {E = K + Vlj + Vcouiomb + 
Vwaiis) for the same configuration. The value of turned 
out to be q^ = 0.238. This last approach mimics the 
nuclear scenario in a more realistic way. 

For both systems we have performed extensive molec- 
ular dynamics simulations following the same approach 
as described in 7]. 

We will first focus on the thermodynamics of these 
systems at a very low densities [p < 0.025cr~'^). For these 
densities a pure LJ system displays a well defined loop in 
its CC as a direct consequence of the ability of the system 
to form well defined fragments in configuration space . 
In the upper panel of Fig.^we show the caloric curves for 
a pure LJ (empty circles), LJ +Coulomb inhomogeneous 
case (full triangles) and homogeneous case (full squares) 
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FIG. 1: Top:Caloric curves for p = O.Olcr"^. Pure 
LJ (empty circles), inhomogeneous (filled squares) 
and homogeneous (filled triangles) systems. Bot- 
tom: Caloric curves for an inhomogeneous system 
with p = 0.015cr-3. A = 147, Z = 62 
(squares) and A — 294, Z = 124 (diamonds). 



for p = O.Olcr^'^. It can be seen that for all the three cases 
a clear loop is displayed. If we restrict our analysis to 
the systems with Coulomb interaction we notice that the 
location of the minima of the loop is shifted with respect 
to the LJ case. This shift is of the order of 0.3e, which 
corresponds to the mean coulombic energy of the system. 
This can be easily verified by plotting the temperature as 
a function oi K + Vlj, in which case the shift dissapears 
denoting that the partition of the total potential energy 
is made into two terms: The LJ potential energy which 
is extremely sensitive to the presence of inner surfaces, 
and the coulomb term, which is, on the contrary, very 
insensitive due to its long range. 

In the lower panel of Fig. we illustrate the fact that 
the presence of the loop is quite resilient to an increase 
of the systems size. In particular we show the CC for 
p — 0.023(T^'^ and A — 294 for the pure LJ system and 
the inhomogeneous case (with Z = 124). It can be seen 
that both the loop and the shift remain present. 

In order to gain further insight into the properties of 
such a system we studied the second moments of the 
distribution of kinetic energies, namely the standard de- 
viation of the kinetic energy per particle and the relative 
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FIG. 2: Top:Relative fluctuations of kinetic energy 
for p — 0.01(7^^. Pure LJ (empty circles), in- 
homogeneous (filled squares) and homogeneous (filled 
triangles) systems. Bottom: Standard deviation 
of kinetic energy per particle for the same cases. 



kinetic energy fluctuations, defined as 



A 



K 



rp2 



(1) 



where N is the number of particles, (t^ the standard 
deviation of the kinetic energy per particle and T the 
temperature of the system. Since kinetic energy fluctua- 
tions and the specific heat are related by jllj 



(2) 



Negative values of the specific heat should be expected 
whenever Ak exceeds the canonical A/f(can) = 1.5 0, 

m 

In Fig. [5] we show these two quantities, Ak (upper 
panel) and (lower panel). Once again the presence 
of the loop in the CC is correlated to the peaks in the 
values of Ak and (Jk- Moreover these peaks are above 
the canonical value (dashed line of Fig. thus denoting 
a negative value of the TRF. It is also interesting to 
notice that for the homogeneous case the size of the peak 
is strongly reduced, which is in accordance with the fact 
that the loop in the corresponding CC is shallower. On 
the other hand for the inhomogeneous case this effect is 
much weaker. 

Following the scheme used in reference |l3j we now ex- 
plore the effect of increasing the density in the above pre- 
sented quantities. The results are summarized in Fig.|31in 
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FIG. 3: Caloric curve for the homogeneous LJ + 
coulomb case for different densities (upper panel). 
Relative fluctuation of the kinetic energy (medium 
panel) (notice that the peak of Ak exceeds the 1.5 
canonical value for p — O.Qla^^). Standard devi- 
ation of the kinetic energy per particle (bottom). 



wich we show the temperature (upper panel), Ak (mid- 
dle panel) and ax (lower panel) for different densities 
(see figure caption for details). In the case of the tem- 
perature it displays a transition from a CC with a loop 
{p = O.OIct"^) to a monotously increasing function for the 
highest density considered (p = 0.60cr~^). For interme- 
diate density values the CC displays a change of slope. 
These features are in complete agreement with the al- 
ready calculated behavior of the CC for the pure LJ case 
|l3j | . The behavior displayed by Ak and aK (middle and 
lower pannels respectively) is also similar to those found 
for the pure LJ case. These results confirm the already 
pointed "triviality" of the inclusion of Coulomb interac- 
tion for small constrained systems. It is worth to mention 
at this point that the behavior of the unnormalized fluc- 
tuation of the kinetic energy aK changes qualitatively 
from displaying a loop to an increasing function of the 
energy at a density between p = 0.3cr~^ and p = OAa~^ 
which has been pointed out jlj| to be the critical density 
for an homogeneous LJ + Coulomb system. 

So far the presence of Coulomb interaction has not 
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FIG. 4: Panels (a) and (6): Local density pro- 
files averages (LDPA) for a liquid-like configuration 
(a) and a gas-like configuration (6). LJ LDPA 
(empty circles), homogeneous LJ + Coulomb LDPA 
(filled triangles), proton LDPA (filled diamonds), and 
neutron LDPA (filled squares). Panels (c) and (d): 
Mass distributions for the LJ (empty circles) and 
homogeneous LJ + Coulomb (filled triangles) cases. 



modified any of the "thermodynamic" properties of the 
system. We now turn to the analysis of its morphological 
properties. We first study the density profiles for the 
three cases studied in this work. For this purpose we 
divide our system in concentrical equally spaced shells 
and calculate the normalized average populations of each 
shell. We have performed this analysis for the low density 
case for two relevant energies, namely E = — 2.0e and 
E — l.Oe for the LJ + Coulomb systems and E = — 2.7e 
and E — 0.7e for the LJ case (panels (a) and {b) in 
Fig. Q . To properly understand this behavior we have 
performed a fragment analysis of the system. We have 
defined fragments according to the following definition 
(see |1,5, 16J): Given a set of particles k, it belongs 

to the same cluster C, if: 



\fieCi, 3jeCi/ Cij < Vm 



(3) 



where aj = V{rij) + (pj — pj)^/Ap, p is the reduced 
mass of the pair {i,j}, and Vmax is an upper limit for 
the potential energy (0 in the LJ case). 

This definition takes into account in an approximate 
way the relative momentum of the particles. The corre- 
sponding mass distributions are displayed in panels (c), 
(d) (for the same energies depicted as before). For the 
sake of clarity we only show the fragment mass spectra 
for the pure LJ and the homogeneous case. We can see 
that for pure LJ and homogeneous case a big fragment 
is formed together with small aggregates giving rise to 
a U-shape distribution so the difference in the density 
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profiles (circles and triangles in panels (a), (5) for LJ, 
LJ + Coulomb respectively) conies from the fact that 
the big fragment for the pure LJ case travels along the 
whole volume remaining mainly in the center of the con- 
tainer whereas for the homogeneous case it remains close 
to the walls due to Coulomb repulsion with the rest of the 
system. For the high energy case the system is mainly 
composed of small fragments, in both analized cases, but 
while for the pure LJ case particles are homogeneously 
distributed inside the volume, for the charged system the 
effect of repulsion gives rise to a higher density close to 
the walls. 

Turning our attention to the inhomogeneous case, it 
is interesting to notice the differences between the pro- 
ton and neutron density profiles (squares and diamonds 
in panels a) and b) of Fig. ^ respectively). For the low 
energy case no major differences are found, which is re- 
lated to the fact that the structure of the big clusters is 
composed by both protons and neutrons. However, the 
scenario changes completely when looking at the system 
at high energies (panel (6)). Since the system is now 
composed mainly of small clusters the density profile of 



the protons is quite similar to the homogeneous case, with 
higher values close to the walls due to Coulomb repulsion, 
but the neutron density profile now clearly resembles the 
LJ profile i.e. an homogeneous distribution inside the 
container. 

To sum up, after analyzing both thermodynamical 
(caloric curves, kinetic energy fluctuations) and morpho- 
logical properties (fragment mass distributions, density 
profiles) of systems with and without Coulomb interac- 
tion we came up with the result that phase transitions 
in small systems, with masses of the order of those rele- 
vant for the nuclear case, are not substantially modified 
by the presence of such long range forces. Nevertheless a 
slight reduction of the signals takes place specially for the 
homogeneous case, which could be useful in the study of 
metallic clusters, where negative specific heats have been 
measured experimentally |l7j . 
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